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We study a Lorentz invariant pairing mechanism that arises when two relativistic spin-1/2 fermions 
are subjected to a Dirac string coupling. In the weak coupling regime, we find remarkable analogies 
between this relativistic bound system and the well known superconducting Cooper pair. As the 
coupling strength is raised, quenched phonons become unfrozen and dynamically contribute to the 
gluing mechanism, which translates into novel features of this relativistic superconducting pair. 
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I. INTRODUCTION 

A large class of superconducting materials can be ac- 
curately described by the BCS theory [l[ , which is based 
upon two major contributions. First, Frolich showed how 
the coupling between electrons and crystal phonons leads 
to an effective attractive interaction between the elec- 
trons 2]. Inspired by this result, Cooper discussed how 
any attractive interaction can bind a couple of electrons 
which lay around a filled Fermi sea 0]. This bound sys- 
tem, known as a Cooper pair, is responsible for several in- 
triguing properties displayed by superconductors, which 
can be described as a many-body coherent state where 
electrons above the Fermi surface are bounded in pairs. 

The oversimplified picture developed by Cooper cap- 
tures the essence of the underlying physical phenomena 
occurring in superconducting solids. In the same spirit, 
we study a simple model of two relativistic fermions with 
an effective attractive interaction. In order to maintain 
the similarities with the Cooper problem, we must fulfill 
the two following requirements: 

Phonon gluing mechanism: In a relativistic scenario, 
the simplest phonon-likc coupling is modeled by a Dirac 
string coupling mechanism where the vibrations of the 
string describe the lattice phonons (see fig. [1] left). This 
interaction is introduced by a non-minimal coupling pro- 
cedure in the free Dirac equation 
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where \ty) stands for the Dirac 4-componcnt spinor, p 
represents the momentum operator, and c the speed of 
light. Here (i — diag(I 2 , — 12), cej = off-diag(<7j, <jj) are 
the Dirac matrices in the standard representation with Uj 
as the usual Pauli matrices Q ■ This Dirac string coupling 
p — > p imw/3r was introduced in 0, 0] as a relativis- 
tic extension of the harmonic oscillator, usually coined 
as the Dirac oscillator, where u> represents the oscilla- 
tor's frequency. In our picture, this frequency effectively 
describes the lattice vibrations and its coupling to the 
fermionic degrees of freedom. 

Two-fermion binding: Regardless of the coupling 
strength, we shall show that such an effective string 
coupling binds relativistic fermions in pairs (see fig. [T] 
right). A Lorentz invariant extension of the Dirac string 



Hamiltonian in Eq. |T]) to two-fermion systems is possible 
@, HI, S El 1 which in the center of mass reference frame 
reads as follows 
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where cti = a ® I4, 0:2 = I4 ® a, (3\ = (3 ® I4, /3 2 = 
I4 (g) (3, and /3i 2 = (3 <E> (3 represent the generalization of 
the Dirac matrices in the two-body Hilbert space. Here 
p := (pj — p 2 )/v2) an d r := (ri — r2)/V2 stand for the 
relative momentum and position operators. 

In this work, we study the binding properties of this 
two-body relativistic Hamiltonian, and discuss under 
which circumstances an analogy to Cooper pairs can be 
performed. Phonons in this relativistic system are dy- 
namical and always provide a pairing mechanism, as we 
shall see. Thus, there is no need to invoke a many- 
body effect through the Pauli principle as in the origi- 
nal Cooper pair scenario. In fact, there are real mate- 
rials which deviate from standard BCS theory. In BCS, 
phonons are quenched and their effect appears as a pair- 
ing energy scale, but they are not explicit in the Hamil- 
tonian. There is an extension of the BCS theory that 
accounts for the effects of dynamical phonons, known as 
the Migdal-Eliashberg theory [H E2, E3| - Our relativistic 
fermionic pairing mechanism is thus closer to this latter 
treatment. 

We shall restrict ourselves to a two-dimensional sys- 
tem, where an exact solution is derived and several in- 
teresting properties can be neatly discussed. Two spa- 





FIG. 1: (left) Electron-phonon interaction in the language of 
Feynman diagrams, (right) Effective electron-electron attrac- 
tive interaction due to the exchange of a lattice phonon. 
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tial dimensions are also natural for other superconduct- 
ing materials like the cuprates [Tjj . In this case, the 
Dirac matrices reduce to the usual Pauli matrices a x = 
<t x , a y = o y , f3 = a z , and the relativistic 1-body state 
can be described by a 2-component spinor. The 2-body 
relativistic Hamiltonian in two dimensions can be written 
as follows 

H 2 u =~^° x ® — I2 ® Vx){Px - imwa z ® cr z a;) 
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II. ENERGY SPECTRUM AND EIGENSTATES 

In two dimensions, chiral creation-annihilation opera- 
tors which carry dual aspects of a left- or right-handed 
symmetry are defined as follows 
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annihilation operators 

„t 



are the usual creation- 
of the harmonic oscillator 



u,^ — (j^Ti — i^wPi ) i an( l A = \J fr/m<-^ is related 
to the ground state width. Using these operators, the 
relativistic Hamiltonian in Eq. ([3]) takes a simpler and 
amenable form 
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stands for the system rest mass, 
2£ is a coupling parameter, and f := 
fiw / trie? controls the strength of the effective interaction. 
Considering the two-body spinorial basis {| tt)>l tl): 
I IT); I II}}: we can understand the Dirac string coupling 
as a four-level system depicted in fig. O 
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FIG. 2: Fermionic spin-flip transitions due to the Dirac string 
coupling and mediated by the creation-annihilation of chiral 
phonons. 

We now proceed to describe the energy spectrum of 
the 2-body interacting relativistic system, in terms of the 



phonon Fock states 
\n r ,ni) := 
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where n r ,ni — 0,1... specify the number of right- and 
left-handed phonons coupling the two-fermion system. 
One immediately observes that the Hilbert space can 
be divided in a series of invariant subspaces TL = 
n =oHn T ni, where each subspace can be described 
by TLn r ni '■= K r .Ti,©K r .ii,- These subspaces are 
spanned by 

ft'n r ,n, =span{|+)|n r ,n/)}, 



=span{| TT)K,«-z + 1), \-}\n r ,m), | ||)|n r + l,ni)}, 
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where the states |-) := (| TT) - I TT))/V2 , and |+) := 
(I TT) + I U*))/V2 are maximally entangled unpolarized 
Bell states. In particular, TL' n n describes a zero-energy 
subspace E+ tTlrini = 0. The Hamiltonian (f5]) in the re- 
maining subspaces TL„ r „, can be expressed as follows 
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where the 2-body interaction couples three different levels 
and can be exactly diagonalized. Using Cardano-Vietta 
solution to third order polynomials, we obtain the follow- 
ing energies 
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These eigenstates are represented for different values of 
the coupling strength £ in Fig. [31 where the chiral quan- 
tum numbers have been set to n r = n; + 1. 

In this figure we observe two different regimes: 
Weak Coupling regime ( < 1: In this case, the low 
energy properties can be accurately described by a two- 
level system. This feature will turn out to be crucial for 
the analogies of the system to a non-relativistic Cooper 
pair discussed in section HTT1 

Strong Coupling regime £ ^> 1 : In this case, the four 
levels become essential in order to describe the low energy 
excitations. Consequently, the description becomes more 
involved but also gives a richer structure that may show 
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certain novel properties with respect to non-rclativistic 
Cooper pairs that are described in section ITVl 

Once the eigenvalues have been obtained, we may de- 
rive the corresponding eigenstates, which we list below 
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Here the indexes j = 1,2,3 correspond to the three differ- 
ent eigenvalues ©. Let us mention that the total angular 
momentum J z := S z + L z is conserved. Thus, the eigen- 
states ([lip have well-defined angular momenta, namely, 
h(n r — tij). Finally, we must consider the consequences 
of fermion indistinguishability. The Symmctrization pos- 
tulate states that a system of identical fcrmions must be 
described in terms of antisymmetrical states, which es- 
tablishes the following constraint 
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(13) 



where P21 stands for the permutation operator that 
swaps the fermion labels 1^2. Considering the eigen- 
states in Eq. pip under the permutation operator, we ob- 
tain the following expressions 
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FIG. 3: Dependence of the two-fermion energy spectrum with 
the dirac string coupling strength. Note that as 1/C, — > we 
approach the strong coupling regime ( 3> 1, whereas l/£ — > 00 
represents a weak attractive coupling (< 1. 



Since these expressions must satisfy the antisymmetric 
condition in Eq. (|13p . the number of chiral quanta are 
constrained as follows 
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Due to the indistinguishability of the relativistic 
fcrmions, the eigenstates \E + t n r ,m) must contain an odd 
number of chiral quanta, whereas \Ej in , rm ) are restricted 
to even number of chiral quanta. 

We have thus derived a complete solution of the rel- 
ativistic Dirac equation for two bodies interacting via a 
Dirac string coupling. Therefore, this 2-fermion system 
belongs to the small class of exactly solvable few-body 
relativistic systems. In sections lllll and ITVl we show that 
this relativistic interaction does indeed lead to the forma- 
tion of bound pairs, both in the weak and strong coupling 
regimes. Furthermore, we present a detailed study of the 
similar properties that the relativistic bound pair shares 
with the well-known non-relativistic Cooper pair. As we 
will see, there are profound analogies in the weak cou- 
pling regime, whereas novel properties are found in the 
strong coupling limit. 



III. WEAK COUPLING REGIME 

The standard description of Cooper pairs in supercon- 
ducting solids is usually performed in a weak coupling 
regime, where a slightly phonon-mediated attractive in- 
teraction binds electron which lay close to the Fermi sur- 
face. We shall consider that the two-body Hamiltonian in 
Eq. effectively describes the gluing mechanism above 
the Fermi sea, and therefore a weak coupling regime is 
obtained when (<C 1. 

In this weak coupling regime, we have seen in Fig. [3] 
that the low-lying excitations can be entirely described 
by a two-level system. This situation is schematically de- 
scribed fig. SI where we see how spin-polarized levels be- 
come decoupled from those responsible of the low-energy 
properties. In this situation, we can obtain an effective 
Hamiltonian for the low energy sector, by adiabatic elim- 
ination. 

Let us consider an arbitrary state ^(f)) € 7Y„ rni 
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whose dynamical evolution, described by the Dirac 
Hamiltonian ([5]), can be represented as 
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In the weak coupling limit, the transitions to the spin- 
polarized {| 11,n r ,ri[ + 1},| U,n r + 1, rij)} upper and 
lower levels can be considered negligible. Therefore, the 
level population does not evolve under the action of the 



two-body interaction = -gjr = 0, and we can adia- 
batically eliminate these two levels. The latter conditions 
substituted in Eq. (fTT)) . give rise to the following relations 
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and an effective two-level dynamics 
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In this sense, we can integrate out the high-frequency 
modes by projecting onto the effective spin-unpolarized 
invariant subspace spanned by H„ n := span{| |j 
, n r , ni), | |t, n r , ti/)}, by means of an orthogonal projec- 
tor 

Kfn, : = I TI,' ? V,n;)(ti.,"-r,n;| + | IT, n r ,ni) (If, n r ,nj | 

(20) 

The z-component of the orbital angular momentum op- 
erator L z = h(a\,a r — ajai) constrained to this invariant 
subspace becomes T"^ ni L z V^ ni = h(n r — ti/)I 2 , which 
allows us to rewrite Eq. (TT9"1) as an effective two-level 
Hamiltonian 
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This effective interaction in the weak coupling regime 
is represented in Fig. 31 where the allowed transitions 
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FIG. 4: Fermionic spin-flip transitions in the weak coupling 
regime. The low energy sector is described by an effective 
two-level system, where spin-flips occur along two different 
channels that include virtual two-phonon transitions and spin- 
polarized states become decoupled. 



can take on two different channels via the consecutive 
creation-annihilation of right- or left-handed phonons. 
This process can be understood as an instance of a su- 
perexchange coupling between the spins | ||) < — > | ||) 
driven by a second order two-phonon process where a 
chiral phonon is virtually created and then annihilated. 
There exist two different exchange paths, as seen in fig. SI 
depending on the left- or right-handed chiralities of the 
virtual phonons involved in the process. 

This effective Hamiltonian (f2Tj) can be exactly diago- 
nalizcd yielding the eigenvalues 
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where the anti-symmetric character of the fermionic 
states has already been considered. Therefore, the low- 
lying solution in the weak coupling regime can be de- 
scribed by the maximally entangled Bell states in the 
spin degree of freedom, and rotational Fock states in the 
orbital degree of freedom. 

Furthermore, these states describe a bound fermion 
pair. In order to show that such binding occurs, we 
must show that the inter-particle distance only attains 
finite values. Let us introduce the square-distance op- 
erator r := x 2 + y 2 , where x :~ (xi — X2) / 'V2 and 
y := (yi — j/2)/v / 2 denote the space coordinate opera- 
tors for the relative fermionic distance. The expectation 
values in the weak-coupling eigenstates (f2"3"| arc 
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which is always finite. We observe the crucial property 
that this system shares with a non-relativistic Cooper 
pair, namely, the pair of relativistic fermions are bounded 
in pairs even for a weak attraction (< 1, 

Another fundamental property that occurs in standard 
Cooper pairs is the presence of an energy gap between the 
paired energy level and the Fermi surface. This energy 
gap is responsible of the stability of Cooper pairs with 
respect to free fermion pairs and is proportional to the 
lattice Debye frequency AE ~ Tvljd- In the relativistic 
regime, we observe that the energy gap with respect to 
the displaced Fermi surface (i.e. e' F = ) is 



Ai? + „ r „ ( — 0, 

AE- nrni = 2hu(n r - 7i/), 



(25) 



and therefore the only stable pair (i.e. AE < 0) is that 
described by the spin-singlet state when 71/ > n r . 

In this sense we obtain a spin-singlet bound pair which 
clearly resembles the situation in standard Cooper pairs 
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where the fermions are also in the singlet state. Further- 
more, we can observe from this discussion that the rela- 
tivistic gap is proportional to the Dirac string frequency 
AE ~ Huj, which plays the role of the usual Debye fre- 
quency in superconducting materials. 

Finally, to take this comparison further, we should 
study the properties of the stable pair eigenstates in 
Eq. (|23p and compare them to the non-relativistic Cooper 
pair features. 

Spin degrees of freedom: In BCS theory, Cooper 
pairs display a singlet state in the spin degree of freedom. 
We observe in Eq. (|23|) that the stable bound fermionic 
pair state has also a spin-singlet component. 

Orbital degrees of freedom: In BCS theory, Cooper 
pairs display a spherically symmetrical wave function 
with an onion-like layered structure. We directly observe 
from fig. [S] that rclativistic bound pair probability distri- 
bution p_ nrTl (r) display a similar spherically symmetric 
onion-like structure. 




-4 -4 



FIG. 5: Spatial probability density profiles for weak coupling 
stable pairs pl*^ (r) := Tv spin ({v\E^ nrni ) {E cS „ rni \r}) with 
7i; > n r . Top figure corresponds to n r = l,n; = 1 proba- 
bility density pf!ii(r), with r = |ri — V2\/V2. Bottom figure 
represents pl^M- 



In this section we have discussed a rclativistic pairing 
mechanism in a weak coupling regime. We have discussed 
in detail several analogies with a non-relativistic Cooper 
pair that naturally arise in this weak coupling limit. Re- 
markably, we obtain binding regardless of the interaction 
strength, which is a fundamental property of BCS sys- 
tems. Additionally, we have shown how the rclativistic 
energy gap scales with the string frequency in the same 
manner as the Cooper pair gap scales with the phonon 
Debye frequency. In this regard, we may conclude that 
the string interaction plays the role of the lattice phonons 
that mediate the effective attraction between fermions in 
the BCS theory. Furthermore, we have also compared 
the nature of the relativistic pair eigenstates with the 
Cooper pair wave functions. We have seen that the rel- 
ativistic bound pair is also described by a spin-singlet 
state and a spherically symmetric onion-like state in the 
orbital degrees of freedom. All these similarities allow 
us to state that this fermionic pairing mechanism can be 
interpreted as a relativistic Cooper pair, since we recover 
most of the usual BCS properties in the weak coupling 
regime. Nonetheless, this Relativistic Cooper pair can 
also be studied in the strong coupling regime, where novel 
properties with respect to the usual Cooper pair in BCS 
theory arise. As we describe below, when the Dirac string 
interaction becomes strong enough, phonons contribute 
dynamically to the gluing mechanism. 

IV. STRONG COUPLING REGIME 

In this section we study the pairing properties of 
the two-body relativistic system in the strong coupling 
regime £ 3> 1. In this limit we must consider the com- 
plete four-level structure of the system ( sec fig. [5]), and 
the energy spectrum becomes clearly more involved in 
Eq. © (see fig. O. 

In fig. [6] we have represented the different energies 
for an interaction strength £ = 5 which lays in the 
strong coupling regime. We clearly see how two levels 
E2,n r ,nt > ^3,n T ,ni become stable pairs with a certain gap 
AE2 : n r , ni < AE 3 , nrtUl < 0. Therefore the strong cou- 
pling gives raise to a couple of stable bound fermionic 
states, namely, 

\E 2 ,n r .n,) ■= -^-[a 2 \-)\n r ,ni) + i/3 2 | Tti*V,»iJ + 1> + 

5 22 

ii, n r + l,rij)]; 
\E 3 ,n r , ni ) := —-[aa\-)\rir,ni) +ifh\ Tt, n r ,n; + 1) + 

"3 

+i<5 3 | U,n r + l,n ; )]. 

(26) 

The spatial probability distribution Pjn r m i r ) of these 
stable fermionic pairs has been represented in fig. [7] in 
the case of n r = ni = 1 . We can clearly observe that the 
density profile preserves the spherically symmetric onion- 
like structure. Nonetheless, noteworthy differences arise 
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with respect to the weak coupling regime ( compare to 
the top fig. [5]). 

Furthermore, these two stable states form a fermionic 
bound pair since the inter-particle distance is finite 

(r> 2 = A 2 [(1 + n r + ni )al + {01 + 5l){2 + n r + n,)] , 

(r) 3 = A 2 [(1 + n r + ni)a\ + (l3 2 3 + S 2 3 ){2 + n r + nj)] ■ 

(27) 

We may conclude that the Dirac string pairing mecha- 
nism leads to bound pairs in the strong coupling regime, 
which display substantial differences with respect to the 
weakly coupled bound states in Eq. (|2U|) . It follows from 
Eqs. (|26|) that the bound pairs are not in a singlet state 
but rather in a linear superposition of different spin sin- 
glet and triplet states entangled with different orbital 
Fock states. In this regard, the relativistic pairing mecha- 
nism does not induce an anti-ferromagnetic ordering any 
longer, and certain spin-polarization may arise depending 
on the value of the coupling strength £. 

It is also instructive to compare the orbital degrees of 
freedom of bound pairs in the weak and strong coupling 
limits. The weakly coupled states in Eq. ((23)) are in or- 
bital Fock states, which represent a certain number of 
vibrational phonons which are frozen in this limit. On 
the other hand, strongly coupled states in Eqs. (|26|) can- 
not be described by a single Fock state, and therefore 
the vibrational phonons acquire a dynamical behavior 
\n r + 1, nj) \n r ,n{) \n r ,ni + 1), which is a clear 
sign of strong coupling in superconductors. We may con- 
clude that the Dirac string phonons, responsible of the 
gluing mechanism, become unfrozen as the coupling bc- 
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FIG. 6: Energy levels of the two-body Dirac oscillator in the 
strong coupling regime £ = 5 as a function of the different 
phonon number. (inset) Detail of the energies corresponding 
to levels E + „ rni , E 3nr7ll 



comes stronger and contribute to the effective attraction 
in a dynamic phenomenon. This is reminiscent of a (s, p)- 
wave symmetry of a SC order parameter. Similar types 
of superconducting states appear in some quantum liq- 
uids like superfluid He 3 : the so-called A- and B-phases 
exhibit different patterns of spin-orbit symmetry break- 
ing IT?| . Layered materials like the ruthenates also 
exhibit unusual symmetry properties like triplet super- 
conductivity [13, 111 [HI . 

We also observe that the strong pairing mechanism 
leads to a couple of possible stable bound pairs (|26|) . 
whereas the weak coupling only produces one stable 
bound pair. Furthermore, the energy gap displayed by 
the bound pairs also depends on the strength of the cou- 
pling. In the weak coupling regime, we have already seen 
that the energy gap scales as AEL ~ hut, whereas the 
scaling in the strongly interaction limit does not present 
such a simple scaling (see fig. [6]). 




FIG. 7: Spatial probability density profiles for strong coupling 
C = 5 stable pairs Pjn r m(r) := Tr apin ({r\Ej nrni ) (Ej nrni \r)) 
with n r = 1, ni = 1. Top figure corresponds to the probability 
density of the stable pair p2ii(r), with r = r-i — r2|/v2- 
Bottom figure represents the probability density of the stable 
pair p3ii(r). 
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V. CONCLUSIONS 

We have studied the relativistic pairing mechanism of 
the two-body Dirac oscillator in two dimensions, where 
a Dirac string coupling leads to fermionic bound pairs. 
We have described two different regimes where binding 
occurs regardless of the interaction strength. 

In a weak coupling regime, the fermionic pair bears 
a strong resemblance to the usual Cooper pair in BCS 
theory. We remarkably found a similar scaling of the en- 
ergy gap, which allows us to identify the Dirac string fre- 
quency to with the lattice phonon Debye frequency lod in 
superconducting materials. Additionally, we found that 
the relativistic bound pair eigenstates are also in a spin 
singlet state, and present a spherically symmetric onion- 
like structure in the probability distribution. All these 
remarkable analogies suggest to interpret this two-body 



Dirac oscillator as an instance of a relativistic Cooper 
pair. Nevertheless, there may also be other types of rel- 
ativistic binding mechanisms yielding also the formation 
of Cooper pairs. 

On the other hand, a strong interaction leads to re- 
markable differences with respect to BCS Cooper pairs. 
In this case, more than one bound pair can be built, 
which in any case is not in a singlet state but rather in 
a linear superposition of singlet and triplet states. Fur- 
thermore, the gluing phonons become unfrozen as the 
coupling strength is raised and dynamically contribute 
to the pairing mechanism. 
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